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The PHENIX experiment at the Relativistic Heavy Ion Collider has measured lu meson produc- 
tion via leptonic and hadronic decay channels in p + p, d+Au, Cu+Cu, and Au+Au collisions at 
y / ijvjv=200 GeV. The invariant transverse momentum spectra measured in different decay modes 
give consistent results. Measurements in the hadronic decay channel in Cu+Cu and Au+Au col- 
lisions show that cj production has a suppression pattern at high transverse momentum, similar 
to that of 7T° and r\ in central collisions, but no suppression is observed in peripheral collisions. 
The nuclear modification factors, 7?aa, are consistent in Cu+Cu and Au+Au collisions at similar 
numbers of participant nucleons. 

PACS numbers: 25.75.Dw 



I. INTRODUCTION 

The measurement of hadrons produced in relativistic 
heavy-ion collisions is a well established tool in the study 
of the hot and dense matter created in the collisions. The 
PHENIX experiment at the Relativistic Heavy Ion Col- 
lider (RHIC) has carried out systematic measurement of 
hadrons in p + p, d+Au, Cu+Cu and Au+Au collisions 
at ^/satjv =200 GeV. When compared to existing mea- 
surements in p+p and d+ Au, measurements in heavy-ion 
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collisions suggest that particle production at high px is 
affected by jet quenching, which is considered to be an 
effect of extremely dense matter created by the collisions 
[l|. High px suppression of 7r° and r/ was measured in 
Cu+Cu and Au+Au @-[| and the nuclear modification 
factors (Raa) of these mesons were found to be consis- 
tent with each other in pt and centrality. A comparison 
with theoretical models was first done for tt° suppres- 
sion in Q , with the result that the suppression increases 
proportional to the number of participating nucleons as 

2 7 3 

N pait . This result is consistent with existing energy loss 
models such as the Parton Quenching Model (PQM) [0]. 

The u> meson comprises light valence quarks similar 
to the 7T° and 77, but has a larger mass (782 MeV) and 
a spin (1). These differences make the omega measure- 
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TABLE I: Summary of the analyzed data samples and uj meson decay channels. Values for previously published PHENIX data 
(PRD83) 0] and (PRC75) @] are given for comparison. 



Data set 



Trigger 



Sampled events 



f Ldt 



Threshold 



Decay channel 



2.74 nb _1 

241 fib' 1 
3.78 pb" 1 

3.06 pb- 1 

813 fib" 1 
80 nb" 1 



Reference 



2003 d+Au 

2004 Au+Au 

2005 p + p 

2005 Cu+Cu 

2007 Au+Au 

2008 d+Au 



ERT 

MB 
ERT 



MB 
ERT 

MB 

ERT 



5.5B 

1.5B 

85B 

8.6B 

5. IB 
160B 



2.4 GeV 

2.4 GeV 

N/A 

0.4 GeV 
1.4 GeV 
1.4 GeV 

N/A 
3.4 GeV 

N/A 

0.6/0.8 GeV 
2.4 GeV 
2.4 GeV 



UJ — > 7T 7T 7T 



UJ — > 7T 7 



UJ — S> 7T 7 



uj — > e T e 
UJ — ?> TT + ir-TV° 
UJ — > 7T°7 



UJ — t 7T 7 
UJ — > 7T (1 7 



UJ — > TV 7 



UJ — S> 7T 7T 7T 
W — > 7T°7 



PRC 75 [8] 

PRC75 [8] 

this work 

PRD83 [7] 
PRD83 [7] 
PRD83 [7] 

this work 
this work 

this work 

this work 
this work 
this work 



merit an additional probe to a systematic study to un- 
derstand mechanisms of parton energy loss and hadron 
production in the collisions. The px dependence of the 
particle production ratio {uj/tt) and the nuclear modifi- 
cation factors (Raa) should add information about the 
parton energy loss mechanism. Furthermore, using mul- 
tiple decay channels: a leptonic channel uj — > e + e~ (with 
branching ratio BR=7.18±0.12x 10 -5 ) and two hadronic 
decay channels w -> n+Tr-n (BR=(89.1±0.7) x 10~ 2 ) 
and ui -> 7r°7 (BR=(8.90+0.27-0.23)xl0- 2 ) @ extends 
the pt range by using the hadronic channels at high px 
and the leptonic channel at low px- 

Baseline measurements of the uj have been performed 
for p+p via the leptonic channel |7], and for the p+p and 
d+Au in the hadronic channel @, [l3] ■ The uj/tt° ratio 
was found to be independent of transverse momentum 
and equal to 0.85 ± 0.05 stat ± 0.09 syst in p + p and 
0.94 ± 0.08 stat ± 0.12 s y st in d+Au collisions for p T > 2 
GeV/c @. 

This article presents the first measurements of uj meson 
production in Cu+Cu and Au+Au collisions at PHENIX 
via the ir°j channel. These measurements permit the 
study of uj suppression at high px- This paper also 
presents measurements of uj in d+Au collisions with sig- 
nificantly reduced uncertainties in the hadronic channel 
and a first measurement in the diclcctronic channel. 



II. EXPERIMENTAL SETUP 

The PHENIX experiment is designed specifically 
to measure electromagnetic probes such as electrons, 
muons. and photons [llj. The detectors of the PHENIX 
experiment can be grouped into three categories: inner 
detectors close to the beam pipe, two central arms with 



pseudorapidity coverage of ±0.35, each covering 90 de- 
grees in azimuthal angle, and two muon detectors, which 
have 2n azimuthal and pseudorapidity coverage of +(1.2- 
2.2) for the south muon arm and -(1.2-2.4) for the north 
muon arm. The central arms are used to measure the uj 
mesons at midrapidity. 

The inner detectors are used for triggering, measure- 
ment of the 2-coordinate of the collision vertex and 
centrality of the interactions with beam-beam counters 
(BBC) and zero degree calorimeters (ZDC). The central 
arms are capable of measuring a variety of particles by 
using Drift Chambers (DC) and Pad Chambers (PC) for 
tracking and momentum measurement of charged par- 
ticles, Ring Imaging Cerenkov detectors (RICH) for the 
separation of electrons up to the tt Cerenkov threshold at 
4 GeV/c, and an Electromagnetic Calorimeter (EMCal) 
for the measurement of spatial positions and energies of 
photons and electrons. The EMCal comprises six sectors 
of Lead Scintillator Calorimeter (PbSc) and two sectors 
of Lead Glass Calorimeter (PbGl). Additional details of 
the PHENIX experimental setup and performance of the 
detector subsystems can be found elsewhere 0, E3] • 

We used data samples collected in 2004, 2005, 2007, 
and 2008 as summarized in Table fl] The data were 
taken using a minimum bias trigger (MB) and the EMC- 
RICH-Trigger (ERT), which is described below. The 
2003 g?+Au data were published in [8| and are included 
here for comparison. The 2005 p + p data were published 
in Q and are used as the baseline for Raa in d+Au, 
Cu+Cu and Au+Au. Two Au+Au data samples were 
taken in 2004 and 2007. The MB trigger required a co- 
incidence between the north and south BBC [llj]. In 
the Au+Au data sample taken in 2004, additional co- 
incidence between the ZDC and BBC was required. To 
enhance the statistics at high px, the ERT trigger was 
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used in p + p, d+Au and Cu+Cu runs which required: 
1) the event to satisfy the MB trigger conditions; 2) the 
presence of at least one high-p^ electron or photon can- 
didate in the event. For electron candidates the ERT 
trigger required a minimum energy deposit of 0.4 (0.6 
and 0.8) GeV/c in a tile of 2 x 2 EMCal towers matched 
to a hit in the RICH in p + p (d+Au) collisions. For the 
photon candidates the ERT trigger required a minimum 
energy deposit of 1.4, 2.4 and 3.4 GeV/c in a tile of 4 x 4 
EMCal towers in p + p, d+Au, and Cu+Cu collisions, 
respectively. In the d+Au and the Cu+Cu analysis, the 
MB data set was used to measure uj production up to 4 
GeV/c in d+Au and 6 GeV/c in Cu+Cu; the ERT sam- 
ple was used at higher pt- The ERT trigger efficiencies 
measured for single photons and electrons and calculated 
for ui mesons is described in Section 1111 Dl 



III. DATA ANALYSIS 

In this section, we describe the event selection and data 
analysis for reconstructing the leptonic (uj — > e + e~) and 
hadronic (uj — > ir + ir~ir° and uj — > tt°j) decay channels of 
the uj. Corrections applied to the raw data to calculate 
the uj meson invariant yields and systematic uncertainties 
related to the measurements are also presented. 



A. Event selection and basic analysis cut 

In the Run-4 PHENIX configuration, the correlation of 
the charge deposited in the BBCs with energy deposited 
in the ZDCs provides a determination of the centrality 
of the collisions. In the other runs, the centralities were 
only determined by using BBC. A Glauber Monte Carlo 
[T4| with the BBC and ZDC responses was used to es- 
timate the number of binary nucleon-nucleon collisions 
(-^Vcoii) and the number of participating collisions (iVp art ) 
for each centrality bin [ll| . 

Events are selected with a reconstructed ^-vertex 
within 30 cm of the center of the interaction region. 
Charged tracks were required to have momenta in the 
range of 0.2 < p T < 5.0 (7.0) GeV/c for the uj -+ e+e" 
analysis in p + p (d+Au) and 0.3 < p T < 8.0 GeV/c 
for the uj — > 7r + 7r~7r° decay channel Q|. Charged parti- 
cles with pt < 0.2 GeV/c have a large bending angle in 
the axial magnetic field of the PHENIX central magnet 
16] and most of them do not pass through the entire 
tracking system. Electrons and positrons are identified 
mainly by the Cerenkov photons emitted in the RICH by 
requiring at least two photomultipliers hit in the RICH 
cells matched to the track [rjj. Also, matching of the 
energy measured for the charged track in the EMCal 
with the momentum measured in the tracking system, 
\E/p— 1| < 0.5, helps to further improve e/n separation. 
Together the RICH and EMCal provide an e/w rejection 
factor of about 1:10 4 . Photon identification is performed 



by the shower shape criteria in the EMCal [llj], and the 
energy of the selected 7 clusters is above 0.2 GeV. 



B. Leptonic analysis 

The leptonic analysis is done only in p + p and d+Au. 
In case of uj — > e + e~, all electrons and positrons recon- 
structed in each event are combined into pairs, resulting 
in signal peaks which sit on top of a combinatorial back- 
ground in the invariant mass distribution. The uncorre- 
cted part of the background is estimated with an event- 
mixing technique, which combines tracks from different 
events with similar event centrality and z-coordinate of 
the collision vertex. Details of the event mixing proce- 
dure are presented in [l9j . 

Figure [1] shows invariant mass spectra of e + e~ pairs 
in p + p and d+Au collisions at ^snn — 200 GeV af- 
ter subtraction of combinatorial background as described 
above. The solid lines show the global fits which in- 
clude: (1) contributions from uj, p and (j) mesons approxi- 
mated with Breit-Wigner functions convolved with Gaus- 
sian distributions to account for the detector mass reso- 
lution; masses and widths of the uj, p and <j) are fixed to 
the PDG values; the p component is calculated assuming 
that ui and p have the same yield and vacuum branching 
ratios; (2) other correlated residual background, which is 
dominated by a contribution from jets, is approximated 
by a second order polynomial function. The detector res- 
olution, which is determined from simulations, is found 
to be dependent on mass and momentum and varies from 
6 MeV/c 2 to 18 MeV/c 2 . 

The u> yield is determined by counting bin contents in 
a 3 ct width (derived from the fitting) and subtracting 
the polynomial background. An associated systematic 
uncertainty from the raw yield extraction is calculated by 
varying the background normalization, fitting functions, 
range and counting methods. The estimated value is 4- 
15% in p + p [7| and 8-15% in d+Au collisions. 

C. Hadronic analysis 

In the uj — » 7T + 7r~7r° and uj — > 7r°7 channels, the first 
analysis step is to reconstruct ir° mesons by combining 
pairs of photons reconstructed in an event. Then the 
mass and width of the ir° peak in the invariant mass dis- 
tribution of photon pairs are parametrized as a function 
of transverse momentum. The 1 er width of the ir° peak 
varies from 13 MeV/c 2 to 9 MeV/c 2 as pr increases from 
1 GeV/c to 4 GeV/c and is determined by the EMCal 
energy resolution. A pair of photons is selected as a 7r° 
candidate if its invariant mass is within 2 a of the recon- 
structed 7T° mass. In Cu+Cu and Au+Au, an additional 
asymmetry cut for ir° candidates is used to reduce combi- 
natorial background, a = \E 7l — E l2 \/\E 11 + E l2 \ < 0.8. 
Selected ir° candidates, which include true 7r°s and com- 
binatorial background are combined either with the third 
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FIG. 1: Invariant mass of e + e~ pairs detected by the PHENIX 
central arms inp + p collisions (left) and minimum-bias d+Au 
collisions (right) at y / sjv"iv=200 GeV and integrated over pr- 
Uncorrelated combinatorial background is subtracted as de- 
scribed in the text. The spectrum is fit to the uj and <f) res- 
onances where the masses and widths are set to the PDG 
values; the Breit-Wigner resonance shape is convolved with 
a Gaussian to account for detector mass resolution estimated 
from simulation and then corrected for the radiative tail. The 
p contribution is shown as the dotted line with an assumption 
that the yield is the same as that of the uj. The residual con- 
tinuum component is estimated by a polynomial fit as shown 
by the dashed line. 



photon with energy E~ f > 1.0 GeV/c for the uj — > 7r°7 
or with a pair of opposite-sign charged tracks for the 
uj — > 7T + 7r~7r° decay. 

In the p+p and d+Au analysis, the uj meson raw yields 
are extracted by fitting the px slices of the invariant mass 
distribution with a combination of a Gaussian for the 
signal and a second order polynomial for the background. 
The width and mass of the reconstructed u> mesons were 
found to be in good agreement with values expected from 
simulation. Details of these analyses are described in [l(| ■ 

In the Cu+Cu and Au+Au analysis, only the uj — > 7r°7 
channel was analyzed due to high combinatorial back- 
ground in the ui —> 7r + 7r~7r° channel. The uncorrelated 
combinatorial background was estimated using an event 
mixing technique where the third photon in the 7r° 7 de- 
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FIG. 2: (a) Invariant mass and scaled mixed background dis- 
tributions for 7r°7 decay at 7 < pr < 12 GeV/c in Au+Au 
collisions (b) Invariant mass distribution after subtraction of 
scaled background. 



cay was taken from the different events with a similar 
centrality and z-vertex. For every px bin the calculated 
background was normalized to match the integral of the 
foreground at an invariant mass Mi nv > 1.75 GeV/c 2 
and then subtracted. An example of the invariant mass 
distribution and normalized background distributions is 
shown in Fig. [^a) with the invariant mass distribution 
after subtraction shown in Fig. HJb). The resulting in- 
variant mass distribution contains residual background 
from correlated particles: the background contributions 
are from K s — > 7r°7r° decays, and ir° and 77, where one 
of the photons from tt° (77) — >■ 77 decay creates a fake 7r° 
candidate for the u> — » 7r°7 reconstruction. The uj — > 7r°7 
peak is further enhanced by a mixed background subtrac- 
tion. Finally, raw yields of uj arc extracted by fitting the 
spectra with a combination of a Gaussian and a polyno- 
mial. The width of the Gaussian is limited ± 1 MeV/c 2 
in the fit to the data as determined from simulation. The 
w yield is calculated as an integral of the Gaussian. 

Systematic uncertainties associated with the raw yield 
extraction are evaluated using different fitting functions 
and ranges, different counting methods and kinematic 
cuts, varying the EMC resolution in simulation, and ap- 
plying different limits for the width of uj peaks in fits 
to data. The estimated value is 13-35% in Cu+Cu and 
20-35% in Au+Au collisions. 



D. Reconstruction efficiencies 

The reconstruction efficiency of the uj is determined us- 
ing a GEANT simulation of the PHENIX detector tuned 
to reproduce the performance of the detector subsys- 
tems. The uj mesons are generated and decayed into cor- 
responding decay channels, and reconstructed with the 
same analysis chain as the real data. The generated uj 
spectra were weighted to match the measured particle 
spectra. It was verified that the simulated positions and 
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FIG. 3: Typical ERT trigger efficiency. Left: trigger efficien- 
cies for single electrons (0.6 GeV threshold) and photons (3.4 
GeV threshold). Right: trigger efficiencies for lo — > e + e~ and 
to — > 7T°7 using corresponding triggered electrons/photons. 



widths of the reconstructed particle peaks are consistent 
with the values measured in real data. 




P T (GeV/c) 



FIG. 4: Typical reconstruction efficiencies for lo — > e + e~ and 
uj — > 7r°7. The curve for ui —¥ 7T°7 includes the embedding 
efficiency in Au+Au collisions: solid, dotted, and dot-dashed 
lines, are respectively for 60-92% , 20-60%, and 0-20% cen- 
trality. 



efficiency e x e cm b X e tr i g , depending on the collision, cen- 
trality, and trigger involved. 

The ERT data sample was used to measure dielectron 
and hadronic decay channels of the u at high px in p + p, 
d+Au and Cu+Cu. The threshold settings for ERT are 
described in Section|TTJ The single particle ERT efficiency 
was measured by dividing the energy spectra of gamma 
clusters or electrons that fired the ERT trigger by the 
energy spectra of all clusters or electrons in the minimum 
bias data sample. A typical example of the ERT trigger 
efficiencies for single electrons and single photons as a 
function of cluster energy is shown on the left in Fig. [3J 
The level of saturation of trigger efficiency curves is below 
100% because of inactive areas of the ERT and the RfCH 
detectors. 

The ERT efficiencies for the uj meson in both the lep- 
tonic and hadronic decay modes were evaluated with the 
help of a Monte-Carlo simulation. For all fully recon- 
structed to mesons, the calculated single photon or elec- 
tron ERT efficiency curves were used to calculate the 
probability that one of the particles in the final state 
fires the ERT trigger. Corresponding trigger efficiencies 
for uj — > e + e~ and uj — > 7r°7 are shown on the right in 
Fig. [3J More detailed descriptions are presented in [7llToj. 

Figure [4] shows typical reconstruction efficiencies e for 
u —> e + e~ and uj — > 7r°7. In the case of Cu+Cu and 
Au+Au collisions, an additional efficiency correction e em b 
due to cluster overlap in high multiplicity environment 
must be applied. In most central Au+Au events, the 
EMCal typically detects more than 300 clusters corre- 
sponding to a detector occupancy of ~10%. To estimate 
the corresponding loss in efficiency, the simulated uj de- 
cays are embedded into real A+A events and analyzed. 
The merging effect results in ^40% loss of reconstruction 
efficiency in 0-20% central Au+Au collisions, ~15% loss 
in 0-20% central Cu+Cu collisions and is almost negligi- 
ble in peripheral collisions. The reconstruction efficien- 
cies derived for Au+Au collisions at different centralities 
are shown in Fig. [5 Finally, in each bin we apply also 
a correction factor [7j to replace the average value of the 
yield in the analyzed px bin by the value of the yield in 
the middle of the bin. 



The reconstruction efficiency is divided into three com- 
ponents: e, etrig, and e e mb- The efficiency e is the re- 
construction efficiency for minimum bias events in a low 
occupancy environment, like in p + p and d+Au colli- 
sions. This efficiency accounts for the limited geometri- 
cal acceptance, resolution and efficiencies of the detector 
subsystems as well as for analysis cuts. When a selec- 
tive ERT trigger is used, an additional trigger efficiency 
factor, etrig, is applied. This factor measures the effi- 
ciency of the ERT trigger logic. For higher multiplic- 
ity collisions, one needs to account for the loss of effi- 
ciency from increased detector occupancy: this is mea- 
sured through the embedding efficiency e em b- A mea- 
sured raw yield then needs to be corrected for the total 



E. Calculation of invariant yields 

In p+p and minimum bias <i+Au collisions, the invari- 
ant yield is related to the invariant cross section as: 



dp* = PP ( dAu) 



I d 2 N 
2ttpt dprdy ' 



(1) 



where <7™ e ' and e^Au are tne total inelastic cross section, 
42.2 mb and 2260 mb respectively. 

For a given centrality bin the invariant yields as a func- 
tion of pt (invariant transverse momentum) are deter- 
mined from: 
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2np T N^ t BR e(p T ) 
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Apr Ay 



(2) 



where N^ t is the number of events for a given centrality 
bin, N(Apx, cent) is the raw yield of oj for each px and 
centrality bin, e(pr), £cmb (pr > cent) and etrigipr) are, as 
previously defined, reconstruction efficiency, embedding 
efficiency and trigger efficiency, respectively. The trigger 
efficiency is applied only for the analyses using the ERT 
data set. BR is the decay branching ratio from Q, (89. 2+ 
OJxKT 2 ) for a; -> tt+tt-tt , (8.90+ 0.27xl0-' 2 ) for a; -> 
7T°7 and (7.16± 0.12xl0~ 5 ) for oj -> e+e - . 



F. Systematic uncertainties 



\/snn—200 GeV are shown in Fig. [5] Measurements were 
performed only in the oj — > 7r°7 channel. Results are pre- 
sented for three centrality bins: 0-20%, 20-60%, 60-92% 
(60-94% in Cu+Cu) and minimum bias collisions. The 
dashed lines represent N co i\ scaled fits to p + p results, 
where N co u values were taken from [IBJ . The results show 
that in peripheral heavy ion collisions oj production gen- 
erally follows binary scaling, while in midcentral and cen- 
tral collisions, production of oj mesons is suppressed at 
high pt- Such behavior is similar to one previously ob- 
served for other light mesons [4|,|20j and can be attributed 
to medium induced effects. 



In addition to uncertainties related to the raw yield ex- 
traction described in the corresponding analysis sections, 
other sources of the uncertainties should also be taken 
into account. Uncertainties of the ERT trigger efficiency 
and acceptance corrections were estimated by varying the 
analysis cuts, energy and momentum scales of the EM- 
Cal and DC by ~1% [7j. Uncertainties of detector re- 
sponse (mainly from the RICH for electron analysis and 
from the EMCal for hadronic analysis) are estimated by 
changing particle identification criteria in the analysis. A 
summary of assigned systematic uncertainties is listed in 
Table [XT] for oj — > e + e~ in p + p and <i+Au and in Table 
Mil for oj — > 7r°7 in Cu+Cu and Au+Au. Those are clas- 
sified into three types: Type A is pr-uncorrelated, Type 
B is ^-correlated and Type C is the overall normaliza- 
tion uncertainty. Total uncertainties for oj — > e + e~ are 
16-24% in p + p and 19-26% in d+Au. The total 
uncertainties for oj — > 7r°7 are 15-37% in Cu+Cu and 
21-37% in Au+Au. Uncertainties for oj — > tt°tt + tt anal- 
ysis are 7-20% in p+p and 10-15% in d+Au, as described 
in 1. 



IV. RESULTS 
A. Invariant transverse momentum spectra 

Figure [5] presents the invariant transverse momentum 
spectra measured for the oj in p+p and d+Au at -^=200 
GeV. Previously published results are shown with open 
markers [sj . Results for different decay channels and data 
samples agree within uncertainties in the overlap region. 
The dashed curves in Fig. [5] are fixed on p + p results 
at px > 2 GeV/c and then scaled by the number of bi- 
nary nucleon-nucleon collisions (N C q\i) estimated using 
Glauber Monte-Carlo simulation 14] for d+Au results. 

Invariant transverse momentum spectra measured for 
the oj meson in Cu+Cu and Au+Au collisions at 



B. 



j /it rati 



Measurement of oj production can be used to study the 
relative production of vector and pseudoscalar mesons 
consisting of the same valence quarks, i.e. oj/tt ratio 
as a function of transverse momentum. In calcula ting 
the oj/tt ratio, the same methodology from 0, [U [22| 
for the 7r + /ir~ and 7r° was used. The charged pion re- 
sults, (ir + + 7r~)/2, were used to extend neutral pion 
measurements at the lower limit of the pt range from 1 
to 0.2 GeV/c. To produce the average pion spectrum in 
p + p [ll| and d+Au collisions [23[ , we simultaneously fit 
(n + +tt~)]2 and ir° spectra with the modified Hagedorn 
function [19j . Inclusion of the charged pion spectrum in 



TABLE II: Summary of assigned systematic uncertainties of 
oj — > e + e~ in p + p and d+Au analysis. 



Source 



p + p 



d+Au 



peak extraction 4-15% (A) 

ERT efficiency 1-3%(B) 

BBC cross section 9.7%(C) 

momentum scale 2-ll%(B) 
acceptance correction 5%(B) 
electron ID 
branching ratio 



8.4-24.1%(A) 
1-7% (B) 
7.9%(C) 
1.2-5.3%(B) 
7%(B) 

10%(B) 
1.796(C) 



TABLE III: Summary of assigned systematic uncertainties of 
oj — > 7r°7 in Cu+Cu and Au+Au analysis. 



Source 



Cu+Cu 



Au+Au 



peak extraction 
ERT efficiency 
energy scale 
energy resolution 
acceptance correction 
conversion 
branching ratio 



13-35% (A) 
3-4%(B) 



4-7% (B) 

2- 3% (B) 

3- 6%(B) 
4.5%(C) 
3.496(C) 



20.1-34.5%(A) 
N/A 
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FIG. 5: Invariant transverse momentum spectra of ui production in p+p and d+Au collisions at -^=200 GeV. The dashed lines 
represent fits to p + p results and those are scaled by the corresponding number of binary collisions for d+Au . The previously 
published PHENIX data (PRD83) @] and tt° (PRC75) @| are shown for comparison. 
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FIG. 6: Invariant transverse momentum spectra of cj production in Cu+Cu (left) and Au+Au (right) collisions from the 
ui — > 7r°7 decay channel for three centrality bins and minimum bias. The dashed lines are the p + p results scaled by the 
corresponding number of binary collisions. The Cu+Cu (left) data was recorded in 2005. For Au+Au (right) the Year 7 and 
Year 4 refer to data taken in 2007 and 2004, respectively. 
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the fit has a small effect in the 1-2 GeV/c overlap region, 
smaller than 5% compared to fitting neutral pions alone. 
The resulting fitted pion distributions are used to calcu- 
late lu/tt ratios for p+p and d+Au. Uncertainties for the 
fit values are evaluated by taking into account statistical 
and systematic uncertainties of the experimental points 
as described in 0, H3] . 

Figure [7] presents the lu/tt ratio measured in p+p colli- 
sions at -\/s=200 GeV as a function of transverse momen- 
tum. Open markers show our previous measurements of 
the lu/tt ratio Q. One can see good agreement between 
previous results and this measurement. For complete- 
ness we also present similar measurements performed in 
lower energy experiments: 7r+Be at ^snn = 31 GeV 
(E706 H), P+P a t %/s=62 GeV (ISR [H). Please note 
that the branching ratio for the lu — > 7r°7 decay was set 
equal to (8.8±0.5)%, which is 6% different from the lat- 
est PDG value of (8.28±0.28)%. Within measurement 
uncertainties the lu/tt ratio in hadronic interactions is 
energy independent at high px- 

A linear fit to the ratio at px > 2 GeV/c gives a value of 
the linear coefficient consistent with zero within less then 
one standard deviation (-0.013 ± 0.009 stat ± 0.014 syst ) 
indicating no significant px dependence of the ratio at 
Pt > 2 GeV/c. A fit to a constant gives a value of the ra- 
tio equal to 0.81 ± 0.02 stat ± 0.09 syst consistent with our 
previous measurement of 0.85 ± 0.05 stat ± 0.09 syst [|. 
The PYTHIA prediction of the ui/tt ratio, shown in Fig. [7] 
with a solid line, lies above the measured ratio. 

The lu/tt ratios measured in minimum bias g?+Au, 
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FIG. 7: Measured lj/tt ratio as a function of px in p + p 
collisions at v / s=200 GeV. (dashed line) Fit of a constant 
value to data points at px > 2 GeV/c. The fit result is 
0.81 ± 0.02 stat ± 0.09 syst . (gray box) The overall er- 
ror of the fitting, (solid line) The pythia prediction [13] 
for p + p at ^=200 GeV. Previously published PHENIX 
results (PRC75) []| and other lower energy experiments at 
y/s^ = 31 GeV (E706) [H and yfs = 62 GeV (ISR) [I| are 
shown for comparison. 



Cu+Cu, and Au+Au collisions at y / ijv r /v=200 GeV are 
presented in Fig. [S] As in the case of p + p colli- 
sions there is no indication that the ratios depend on 
transverse momentum for pt > 2 GeV/c. Fits to 
a constant for pt > 2 GeV/c give the following val- 
ues of the lu/tt ratio: 0.75 ± 0.01 stat ± 0.08 syst in 
d+Axi, 0.71 ± 0.07 stat ± 0.07 syst in Cu+Cu and 
0.83 ± 0.09 stat ± 0.06 syst in MB Au+Au collisions. 
Within the uncertainties the lu/tt ratios measured in dif- 
ferent collision systems for px > 2 GeV/c are in agree- 
ment. This agrees with previous measurements in d+Au 
Q within the uncertainties. The ratios in various col- 
lision systems imply similar suppression factors and px 
dependences within the uncertainties for the lu and tt 
production in nucleus- nucleus collisions at high px- 



C. Nuclear modification factors 

To quantify medium induced effects on high px particle 
production the nuclear modification factor is defined as: 



Rab(pt) = 



d 2 N AB /dydp T 



((N coll )/aZ el ) 



d 2 a pp /dydpx ' 



(3) 



where d 2 Np^/dydpx is the differential yield per event in 
nucleus-nucleus collisions, (iVcoii) is the number of binary 
nucleon-nucleon collisions averaged over the impact pa- 
rameter range of the corresponding centrality bin calcu- 
lated by Glauber Monte-Carlo simulation 14i] , and a™ el 
and d 2 <7p P /dydpx are the total and differential cross sec- 
tions for inelastic p + p collisions, respectively. In the 
absence of medium-induced effects, the yield of high-p T 
particles is expected to scale with (N co n), resulting in 
i?AB=l at high-pT- 

Figure [5] presents RdAu measured for the lu in mini- 
mum bias, most central and peripheral d+Au collisions at 
y/sNN=200 GeV. Good agreement is observed between 
different decay modes, and between new and previously 
published PHENIX lu results Q shown with open mark- 
ers. For comparison we also present tt° results published 
in [231 ] in the figure. In peripheral collisions the measured 
values of RdAu are consistent with unity over the whole 
px range of measurements. In most central collisions a 
modest Cronin-like enhancement is observed in a range 
of px from 2 GeV/c to 6 GeV/c and suppression of lu 
production at 

px > 8 GeV/c. A similar enhancement at 2-6 GeV/c 
was previously observed for neutral and charged pi- 
and 



28j . Suppression of lu pro - 



mcsons 



ons [21], |23 

duction at higher px is in agreement with tt° results [23jj . 
Similarity of the observed effects for the mesons with 
very different masses suggests that the collective nuclear 
effects occur at the partonic level (29l - [3l| . 

Nuclear modification factors measured in Cu+Cu and 
Au+Au collisions at ^/snn— 200 GeV as a function of 
px are shown in Fig. 1101 Results are presented for min- 
imum bias, most central (0-20%), midcentral (20-60%) 
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(0-88%) for uj -> e+e", tt'V 

versus transverse momentum in Cu+Cu (0-94%) for ui 
and in Au+Au (0-92%) for uj — > n ^. The dashed lines and 
boxes are a fit of a constant value to the data points at px > 
2 GeV/c in p + p (Fit result: 0.81 ± 0.02 stat ± 0.09 syst ). 
The previously published data (PRC75) [ij are shown for 
comparison. 



FIG. 9: Nuclear modification factor, RdAu, measured for the 
uj in 0-88, 0-20, and 60-88% centrality bins in d+Au col- 
lisions at y / s=200 GeV. The box at the right edge of the 
constant fit line shows the uncertainty of the fit. The previ- 
ously published data for uj (PRC75) l| and tv° (PRL98) [H 
are shown for comparison. 
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FIG. 10: Raa of the uj in Cu+Cu (left) and Au+Au (right) collisions from the uj — > it°j decay channel for three 
centrality bins and minimum bias. The uncertainty in the determination of p + p scaling is shown as a box on the 
left in each plot. Rhombuses in each plot are 7?aa of n° in Cu+Cu [22] ] and Au+Au [|| shown as a comparison. 
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FIG. 11: 7?aa for the uj meson integrated over the range pr > 
7 GeV /c as a function of the number participating nucleons 
(AW). Results for tt° (PRL101) @] and n (PRC82) [| are 
shown for comparison. The dashed line shows the fitted frac- 
tional energy loss function, Raa = (1 — SoNp a , Tt ) n ~ 2 . 



and peripheral (60-94% in Cu+Cu; 60-92% in Au+Au) 
collisions. The nuclear modification factors do not de- 
pend on pt for pt > 6 GeV/c at all centralities. For 
Np aT t > 34 suppression of uj production begins to be ob- 
served, with suppression increasing as N part increases. 

Figure [11] shows Raa values integrated for px > 7 
GeV/c as a function the number of participants. For 
u) mesons we present four centrality bins in <i+Au, and 
three centrality bins in Cu+Cu and Au+Au. For compar- 
ison the average values of Raa for n° .4] and r\ mesons 
[j| mesons for pt > 7 GeV/c are also plotted. To see 
whether the u> follows the suppression pattern of 7r° and 
77, the integrated Raa vs iV par t dependence is fit to a 
fractional energy loss function Raa = (1 — 5oAp art ) n-2 
0, [Uj]. The parameter n, which is an exponent of the 
power law fit to the uj pr spectrum measured in p + p for 
Pt >5 GeV/c 0, was fixed to 8. The fitting gives x 2 /ndf 
less than three and parameters So = (9.9+0.7) x 10 -3 and 
a = 0.55 ± 0.01. As in [4( we find the parameter a con- 
sistent with predictions of the GLV 33] and PQM @ 
models (a ~ 2/3). Therefore, we can conclude that uj 
production has a similar suppression pattern as ir° and 
r\ which supports the scenario that the energy loss takes 
place at the parton level in the hot and dense medium 
formed in the collisions. 



ate pr 2-6 GeV/c and suppression for pt > 8 GeV/c in 
most central d+Au collisions. The measurement of the 
nuclear modification factor for the w meson in Cu+Cu 
and Au+Au collisions show that uj production has a sim- 
ilar suppression pattern as the 7r° and 77 within model 
agreement, thus supporting the scenario that the energy 
loss takes place at the partonic level in the hot and dense 
medium formed in the collisions. 
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V. SUMMARY 



We measured production of the ui meson via both 
leptonic and hadronic decay channels in p + p, <i+Au, 
Cu+Cu and Au+Au at ^/sjvat=200 GeV. The invariant 
transverse momentum spectra show good agreement in 
different decay channels in p + p and d+Au. The RdAu 
shows a moderate Cronin like enhancement at intermedi- 
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